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ABSTRACT 

Recent stellar evolutionary calculations of low-metallicity massive fast-rotating main-sequence stars 
yield iron cores at collapse endowed with high angular momentum. It is thought that high angu- 
lar momentum and black hole formation are critical ingredients of the coUapsar model of long-soft 
7-ray bursts (GRBs). Here, we present 2D multi-group, flux- limited-diffusion MHD simulations of 
the collapse, bounce, and immediate post-bounce phases of a 35-Mq collapsar-candidate model of 
Woosley & Heger. We find that, provided the magneto-rotational instability (MRI) operates in the 
differentially-rotating surface layers of the millisecond-period neutron star, a magnetically-driven ex- 
plosion ensues during the proto-neutron star phase, in the form of a baryon-loaded non-relativistic 
jet, and that a black hole, central to the coUapsar model, does not form. Paradoxically, and although 
much uncertainty surrounds stellar mass loss, angular momentum transport, magnetic fields, and the 
MRI, current models of chemically homogeneous evolution at low metallicity yield massive stars with 
iron cores that may have too much angular momentum to avoid a magnetically-driven, hypernova-like, 
explosion in the immediate post-bounce phase. We surmise that fast rotation in the iron core may 
inhibit, rather than enable, coUapsar formation, which requires a large angular momentum not in the 
core but above it. Variations in the angular momentum distribution of massive stars at core collapse 
might explain both the diversity of Type Ic supernovae/hypernovae and their possible association 
with a GRB. A corollary might be that, rather than the progenitor mass, the angular momentum 
distribution, through its effect on magnetic field amplification, distinguishes these outcomes. 
Subject headings: MHD - stars: neutron - stars: supernovae: general - neutrinos - rotation - Gamma- 
ray: bursts 



1. INTRODUCTION 

There is mounting observational evidence for the as- 
sociation between long-soft 7-ray Bursts (GRBs) and 
broad-lined Type Ic supernovae (SNe; see Woosley & 
Bloom 2006 for a review). Such hydrogen-deficient (and, 
perhaps, also helium-deficient) progenitors are compact 
and, if fast rotating in their core at collapse, fulfill critical 
requirements for the formation of a coUapsar (Woosley 
1993). The engine that converts energy from long-term 
accretion of disk material onto the black-hole (BH) may 
power a relativistic jet in the excavated polar regions. 
The jet breaks out of the progenitor surface while equa- 
torial accretion continues. Depending on the BH mass 
and the angular momentum budget in the collapsing en- 
velope, this "engine" may operate for seconds, i.e. as long 
as typical long-soft GRBs. Accompanying this beamed 
relativistic polar jet might be a disk wind, fueled by neu- 
trinos or MHD processes, which would explode the Wolf- 
Rayet envelope. This explosion and the radioactive ^^Ni 
material produced might lead to very energetic, broad- 
lined, Type Ic SN of the hypernova variety (Iwamoto 
et al. 1998; MacFadyen & Woosley 1999; Hjorth et al. 
2003; Stanek et al. 2003). 

State-of-the-art radiation-hydrodynamics simulations 
including a sophisticated equation of state (EOS) and 
detailed neutrino transport (Buras et al. 2003; Bur- 
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rows et al. 2006,2007a; Kitaura et al. 2006; Marek 
& Janka 2007; Mezzacappa et al. 2007) suggest that 
while the neutrino mechanism of supernova explosions 
may work for the lower-mass massive progenitors, it may 
not for the more massive progenitors, characterized by an 
ever higher post-bounce accretion rate onto the proto- 
neutron star (PNS). Burrows et al. (2006,2007a) have 
suggested that an acoustic mechanism will work for all 
slowly rotating progenitors that do not explode by other 
means within the first second after bounce. However, 
massive star cores endowed with a large angular mo- 
mentum at the time of collapse should experience the 
magneto-rotational instability (MRI; Balbus & Hawley 
1991; Akiyama et al. 2003; Pessah et al. 2006; Shibata 
et al. 2006; Etienne et al. 2006), with the potential 
to exponentially amplify weak initial fields on a rota- 
tion timescale. The saturation values of such fields are 
ultimately set by the free-energy of differential rotation 
available in the surface layers of the PNS (Ott et al. 
2006), and can be large, i.e., on the order of 10^^ G at 
a radius of a few tens of kilometers. The corresponding 
magnetic stresses at the neutron star surface lead sys- 
tematically to powerful jet-like explosions ~100ms after 
bounce (see, e.g., Ardeljan et al. 2005; Yamada & Sawai 
2004; Kotake et al. 2004; Sawai et al. 2005; Moiseenko 
et al. 2006; Obergaulinger et al. 2006; Burrows et al. 
2007b, hereafter B07; Dessart et al. 2007). 

In this letter, we investigate, in the context of the col- 
lapsar model, the potential implications of this magnetic 
explosion mechanism. Our study focuses on the imme- 
diate post-bounce phase, whose importance was empha- 
sized by Wheeler et al. (2000,2002). This is in contrast to 
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TABLE 1 

Properties of our two MHD-VULCAN/2D simulations of 
THE 350C collapsar model of WH06. 
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Note. — icnd gives the time at the end of each simulation, while to 
is the time when the rate of polar mass ejeetion first overcomes equa- 
torial mass accretion. All quoted quantities in the table correspond to 
the final time in each simulation, while times are given with respect to 
core bounce. Mio (Pio) corresponds to the total baryonic mass (av- 
erage rotation period) inside the lO^^gcm"^ isodensity contour. Egas 
(Eg g) is the Bernoulli (Foynting) power in the ejecta, obtained by 
integrating the corresponding flux over a shell with a radius of 500 km. 
[Sec text for additional information.] 

previous work which explored only the phase subsequent 
to BH formation (MacFadyen & Woosley 1999; Aloy et 
al. 2001; Zhang et al. 2003; Proga 2005). Indeed, two 
terms sometimes used in the collapsar context are "failed 
SN" (MacFadyen & Woosley 1999) and "prompt BH for- 
mation" (MacFadyen et al. 2001). Our analysis supports 
the idea that the conditions for the collapsar model, as 
stated so far, are also suitable for a magnetically-driven 
explosion in the immediate post-core-bounce PNS phase, 
and that BH formation may be so delayed for a range of 
putative progenitor models that it does not in fact oc- 
cur"^. In §2, we present radiation MHD simulations with 
the code VULCAN/2D (Livne et al. 2004,2007) of a 
collapsar-candidate model that support this thesis. In 
§3, we discuss the implications of our results for stellar 
evolutionary models that might lead to coUapsars and/or 
hypernovae. 

2. MODEL AND RESULTS 

We present results from two-dimensional, rotat- 
ing, multi-group, flux-limited diffusion magneto- 
hydrodynamics simulations, using VULCAN/2D (Livne 
et al. 2004,2007; see also appendices in Dessart et al. 
2006 and Burrows et al. 2007a), of the 350C progenitor 
model and collapsar candidate of Woosley & Heger 
(2006; hereafter WH06). The numerical procedure we 
follow is identical to that of B07 in every respect, except 
the choice of progenitor. WH06's model is evolved from 
a 35 Mq zero-age main sequence star endowed with a 
total angular momentum of 1.4x 10^'^ erg • s, a metallicity 
of 1% the solar value, and a reduction by a factor of 10 
in the prescribed mass loss rates during the Wolf-Rayet 
phase. Our simulations extend out to a maximum 
radius of 5000 km (which contains ~3Mq) and cover 
a 90° quadrant, bounded by the rotation axis and the 
equator. We adopt WH06's initial rotational, density, 
temperature, and electron-fraction profiles for that 
progenitor. For the magnetic field distribution, we start 
with magnitudes and morphology that are consistent 
with the 350C model of WH06. We use a field uniform 

^ In the present context, BH formation is never prompt, since it 
takes a finite time, on the order of seconds, for the PNS to accu- 
mulate the critical mass at which it experiences the gravitational 
instability. This is in contrast with super-massive stars, such as 
the progenitors of pair-instability SNe, which may form an appar- 
ent horizon during collapse and thus "directly" transition to a BH 
(Liu et al. 2007). 



within 3000 km, and dipolar beyond. In our reference 
model, MO, we employ initial poloidal and toroidal field 
magnitudes of 2x10^° G and 8x10^^ G, respectively, in 
close quantitative agreement with the 350C model of 
WH06. However, we also study a model. Ml, with an 
initial poloidal field that is five times stronger. This 
leads to a magnetic field energy at ~ 100 ms after core 
bounce that is closer to the value expected at the PNS 
surface, were we to adequately resolve the MRI (B07). In 
Table 1, we give important quantities characterizing the 
two simulations performed. Note that if magnetic fields 
are ignored in the pre-collapse evolutionary calculations 
of WH06, core angular velocities reach 5-22rads~^, 
much larger than the 1.98rads~^ achieved in the 350C 
model. A bounce at sub-nuclear densities may ensue 
and lead to BH formation (Akiyama & Wheeler 2005). 
Ignoring magnetic torques in the models most prone to 
magnetic-field generation during the pre-collapse phase 
seems inconsistent. Thus, we focus on the more slowly 
rotating progenitors, evolved with magnetic fields, which 
inevitably bounce at nuclear densities. 

After an initial collapse phase that lasts '^245 ms, the 
central density reaches ~3xl0^'^gcm~^, the EOS stiff- 
ens, a shock is born, and propagates outward, but is de- 
bilitated by the photo-dissociation of the infalling outer 
iron-core and the burst of electron neutrinos. The shock 
stalls at ^^150 km, and within a few tens of ms after 
bounce, it becomes increasingly aspherical. The net gain 
from neutrino emission and absorption processes, the en- 
tropy, and the material accretion rates get progressively 
larger at larger latitudes as the degree of oblateness of the 
fast-rotating PNS increases. Subsequent to the amplifi- 
cation due to compression by a factor of ~2500 in both 
magnetic field components, the toroidal magnetic field 
increases after bounce due to the winding of the poloidal 
field component. At the same time, accretion of the outer 
magnetized core continues and enhances the total mag- 
netic field energy. By 150 ms (300 ms) after bounce, the 
magnetic pressure at the surface of the PNS along the 
pole in the Ml (MO) model is comparable to the gas pres- 
sure, and a bipolar, magnetically-driven, baryon-loaded, 
and non-relativistic jet is initiated, reversing accretion 
into ejection along the polar direction. As shown in the 
left panel of Fig. [U the initial jet mass-loss rate is only 
^O.OIMq s~^ but in the Ml model and by ~350ms af- 
ter bounce it exceeds the accretion rate. At this time, the 
accumulated baryonic PNS mass is only 1.93 Mq (middle 
panel of Fig.[T|), and, thus, well below the 2.17M0 bary- 
onic mass transition to BH formation that we derive with 
our Shen EOS (Shen et al. 1998). Note that for a solid- 
body rotator this limit may increase by up to 10-20% 
(Cook et al. 1994), but for a differentially-rotating neu- 
tron star this limit may be considerably larger (Baum- 
garte et al. 2000). Hence, rotation, and in particular 
differential rotation, enhances the potential for explosion 
during a PNS phase. By 650 ms after bounce, the explo- 
sion energy in the Ml model reaches ^3.3 B (10^^ erg = 
1 Bethe [1 B]) (see right panel of Fig.[T]), although due to 
the continued accretion along near-equatorial latitudes 
a quasi-steady-state is reached with an explosion power 
sustained at ~10Bs~^. As shown in Fig. [21 the jet re- 
sembles a magnetic tower (Lynden-Bell 2003; Uzdenski 
& MacFadyen 2001), but is confined primarily by the 
ram pressure of the infalling dense envelope (B07). As 
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Time afler bounce [s] Time after bounce [s] Time after bounce [s] 



Fig. 1. — Left: Time evolution of the instantaneous integrated mass flux accreting {solid line) and outflowing {dashed lines) through a 
shefl at a radius of 500km, for models MO (black) and Ml (red). Middle: Same as left, but for the total mass interior to the iso-density 
contour corresponding to lO^'^gcm"^. Right: Same as left, but for the explosion energy {solid line) and the net integrated neutrino gain 
{dashed lines) outside the high-density regions bounded by the lO^^gcm"^ contour. [See text for discussion.] 



time progresses, its base broadens and, given the quasi- 
steady jet conditions, the mass ejection rate grows and 
the accretion is hmited to progressively smaller latitudes. 
Both cause and effect, the increase in ejecta volume en- 
hances the neutrino contribution to the explosion energy, 
although it remains a subdominant part of the total by 
the end of the simulation. Extraction of core rotational 
energy by magnetic torques is also in evidence in the Ml 
model from the increase in the average PNS rotation pe- 
riod^ from 8 to 12 ms between 200 ms and 600 ms after 
bounce, while over half this interval the average period 
decreases from 5 to 4 ms in the weakly exploding MO 
model. The decrease in the free energy of rotation in 
the Ml model is on the order of 3B, and is comparable 
to the magnitude of the explosion energy. This supports 
the idea that core rotation energy fuels the magnetically- 
driven ejecta. Hence, the Ml model, modified slightly to 
yield fields at saturation that agree roughly with what 
would obtain in the presence of the MRI, boasts a clear 
and powerful explosion. In this model, and once the ex- 
plosion is well established, the PNS loses mass at a steady 
rate, and has a mass of only ~1.7M0 at the end of the 
simulation. The broadening of the base of the jet sug- 
gests that the explosion will not choke (nor induce any 
significant fallback) , and encompassing a larger solid an- 
gle, will instead lead to explosion in all directions. Hence, 
such an object is unlikely ever to transition to a BH and 
to lead to a collapsar. 

By contrast, in the MO model, the explosion emerges 
later, when the neutron star baryonic mass has already 
accumulated ~2.1M0, and thence may be susceptible 
to collapse to a BH. The free-energy of core rotation has 
been partially tapped, but the potential subsequent pow- 
ering of a GRB may not be compromised. This model, 
by mimicking more slowly rotating cores or an ineffi- 
cient MRI, offers a limiting case for the formation or 
non-formation of a BH, and a possible collapsar. 

3. DISCUSSION 

The potential for exponential growth on a rotational 
timescale of initial seed magnetic fields by the MRI (Shi- 
bata et al. 2006; Etienne et al. 2006), fueled by the free 
energy of core rotation, makes the initial angular mo- 
mentum budget of the progenitor star the key parameter 
in determining the outcome during the immediate post- 

* We define the average rotation-period as the period of the 
rigidly-rotating PNS that has the same total angular momentum 
and structure inside the 10^'^ gcm~'^ isodensity contour. 



bounce phase (B07). A magnetically-driven, baryon- 
loaded, and non-relativistic, explosion is obtained for 
WH06's 350C collapsar candidate model, evolved at low 
metallicity from a 35 Mq fast-rotating main-sequence 
star. The explosion occurs ~200ms after bounce and 
reaches ^3B ~400ms later. After an initial accretion 
phase, the steadily decreasing PNS mass reaches only 
~1.7M0 at the end of the simulation, and, thus, the 
quasi-steady explosion we observe suggests that BH for- 
mation is unlikely to occur. Moreover, baryon contam- 
ination prevents the ejecta from becoming relativistic. 
Note that the production of a GRB in the collapsar con- 
text is contingent on the gravitational collapse of the 
PNS to a BH. 

The recent stellar evolutionary calculations of Yoon 
& Langer (2005), WH06, and Meynet & Maeder (2007) 
of fast-rotating main-sequence objects at low metallic- 
ity systematically predict such fast-rotating cores at col- 
lapse. Starting from similar conditions for a 35-40 M© 
star, but using different mass-loss "recipes," they obtain 
very similar rotational profiles in the inner core. Allow- 
ing for anisotropic mass loss (Meynet & Maeder 2007), 
a model of C. Georgy (2007, priv. comm.) suggests an 
even larger (by a factor of two) specific angular momen- 
tum in the inner 3 M© at the end of silicon core burning. 
Despite the agreement between these different evolution- 
ary computations, the magnetically-driven explosion and 
the "failed" BH formation described here are conditional 
on the uncertain treatment of mass loss, angular mo- 
mentum transport, and magnetic processes (Spruit 2002) 
during the pre-coUapse evolution. 

At very low metallicities, radiatively-driven winds of 
massive stars are inhibited by the lack of metals (Ku- 
dritzki 2002; Vink et al. 2001; Vink & de Koter 2005), 
whose optically-thick lines intercept radiation momen- 
tum (Castor et al. 1975). Recent revisions down- 
ward of mass-loss rates due to clumping (Owocki et al. 
1988; Bouret et al. 2005; FuUerton et al. 2006) sug- 
gest, however, the potential importance of episodic out- 
bursts, akin to the 1843 giant eruption of Eta Carina 
(Smith & Owocki 2006). The metallicity dependence of 
such phenomena is entirely unknown, mostly because the 
fundamental cause of the outburst remains a mystery. 
While line driving seems excluded, continuum driving of 
a porous medium at super-Eddington luminosities has 
been proposed by Owocki et al. (2004) as an alterna- 
tive. Finally, mass loss in fast-rotating, and sometimes 
critically-rotating (Townsend et al. 2004), envelopes 
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Fig. 2. — Colormap of the entropy at 666 ms after bounce 
for model Ml, overplotted with white iso-density contours (every 
decade downward from lO^'^gcm"^) and velocity vectors (length 
saturated to 15% of the width of the display and corresponding to 
a velocity of 30,000 km s"!.) 

is complicated by the effects associated with centrifu- 
gal support, surface oblateness, and gravity-darkening 
(Cranmer & Owocki 1995; Owocki et al. 1996), so that 
the mass-loss "recipes" used in stellar evolutionary mod- 
els are not always substantiated by observational and 
theoretical evidence. At present, and in light of our sim- 
ulations, it appears that chemically-homogeneous evolu- 
tion of fast-rotating main-sequence massive stars at low 
metallicity systematically yields iron cores at collapse 
that may have too much angular momentum, a property 
that prevents the formation of a coUapsar. Uncertainties 
in the modeling of the pre-coUapse evolution may result, 
however, in slower-rotating iron cores^ and, thus, might 
inhibit an early magnetically-driven explosion in favor of 
black hole, and perhaps coUapsar, formation. 

We conclude that variations in the angular momentum 
distribution of pre-coUapse massive stars may lead to dif- 
ferent post-bounce scenarios. Non- or slowly-rotating 
progenitors may explode with weak/moderate energy 
(<1 B) through a neutrino or an acoustic mechanism <1 s 
after bounce, or may collapse to a BH. Objects with large 
angular momentum in the envelope, but little in the core, 
may proceed through the PNS phase, transition to a BH 
and form a coUapsar with a GRB signature. Owing to 
the modest magnetic- field amplification above the PNS, 
a weak precursor polar jet may be launched, soon over- 
taken by a baryon-free, coUimated relativistic jet. At the 
same time, the progenitor envelope is exploded by a disk 

^ Note that the rotational energy i?rot is a stiff function of an- 
gular velocity w, i.e., E^ot <x ui^. 



wind, resulting in a hypernova-like SN with a large lu- 
minosity (large ^^Ni mass). Finally, and this is what we 
conclude here, objects with large angular momentum in 
the core may not transition to a BH. Instead, and fueled 
by core-rotation energy, a magnetically-driven baryon- 
loaded non-relativistic jet is obtained without any GRB 
signature. The explosion has the potential of reach- 
ing energies of a few B to 10 B, and for viewers along 
the poles of looking like a Type Ic hypernova-like SN 
with broad lines. For a viewer at lower latitudes, the 
delayed and less energetic explosion nearer the equator 
may look more like a standard Type Ic SN (Hoflich et 
al. 1999). This volume- restricted jet-like explosion is 
dimmer, as the amount of processed ^^Ni may be signif- 
icantly less than the ^0.5 obtained in the coUapsar 
context (MacFadyen & Woosley 1999). Hence, magnetic 
processes during the post-bounce phase of fast-rotating 
iron cores offer a potential alternative to coUapsar for- 
mation and long-soft GRBs by producing non-relativistic 
non-Poynting-flux-dominated baryon-loaded hypernova- 
like explosions without any GRB signature. Importantly, 
while our study narrows the range over which the coUap- 
sar model may exist, it also offers additional routes to 
explain the existence of GRB/SN-hypernova events like 
SN 1998bw (Woosley et al. 1999), and hypernova events 
like SN 2002ap without a GRB signature (Mazzali et al. 
2002). 

More generally, magnetic effects should naturally arise 
in the context of gravitational collapse and fast rotation. 
The resulting angular momentum of newly-formed BHs 
and magnetars, for example, would be reduced, perhaps 
considerably, by any prior magnetically-driven explosion, 
and, thus, may decrease the power of subsequent mass 
ejections from compact objects (see, e.g., Thompson et 
al 2004). 
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